% main_fig1.m
clear; clc; close all;

% 参数设置
mu_values = linspace(0, 4, 1000);  % μ 参数从 0 到 4，共 1000 个点
iterations = 1000;                 % 总迭代次数
last = 100;                      % 仅绘制最后 100 次迭代的结果

hold on;
for mu = mu_values
    x = logistic(mu, 0.5, iterations);
    % 绘制当前 μ 对应的最后几次迭代结果
    plot(mu * ones(1, last), x(end-last+1:end), '.k', 'MarkerSize', 1);
end
xlabel('\mu');
ylabel('x');
title('Logistic 映射分岔图');
grid on;
hold off;
